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ABSTRACT 

A long XMM-Newton exposure is used to observe solar wind charge exchange (SWCX) emis- 
sion from exospheric material in and outside Earth’s magnetosheath. The light curve of the O VII 
(0.5 — 0.62 keV) band is compared with a model for the expected emission, and while the emission 
is faint and the light curve has considerable scatter, the correlation is significant to better than 
99.9%. This result demonstrates the validity of the geocoronal SWCX emission model for predict- 
ing a contribution to astrophysical observations to a scale factor of order unity (1.36). The results 
also demonstrate the potential utility of using X-ray observations to study global phenomena of 
the magnetosheath which currently are only investigated using in situ measurements. 

Subject headings: X-rays: general, X-rays: ISM, solar system: general, interplanetary medium, (Sun:) 
solar wind 

striping. However, it was during the ROSAT All- 
Sky Survey (RASS, Snowden et al. 1995, 1997) 
that the insidious nature of the solar system emis- 
sion became clear, even if the source was not yet 
identified. The RASS data had considerably bet- 
ter statistics and coverage than previous surveys 
allowing the striping to be much more obvious. 
The contamination is clearly visible in the un- 
cleaned maps of the j keV band (Figure 1, upper 
panel). At the time the source of the striping was 
unknown and was treated empirically by modeling 
it as a temporally-varying flat field and then sub- 
tracting it from the data (Figure 1, lower panel). 
In practice this worked relatively well although 
some residual striping remained. For lack of a 
better identification and because of the typical 
relative durations (e.g., compared to auroral X- 
rays which had durations of order 10 minutes), 
the count-rate excesses responsible for the strip- 
ing were referred to as “long-term enhancements” 
(LTEs, Snowden et al. 1994, 1995). 

ROSAT pointed observations were also often af- 
fected by this contamination and methods were 
developed to subtract them from images (Snow- 


1. Introduction 

Diffuse X-ray emission from the solar sys- 
tem has been observed for several decades as an 
unidentified contamination component in observa- 
tions of the soft X-ray background. It was seen in 
the Wisconsin sounding-rocket survey data as off- 
sets between overlapping parts of adjacent fields 
(McCammon et al. 1983; Burrows 1982), and in 
the HEAO-1 A2 maps (Garmire et al. 1992) as 
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den et al. 1994). One such observation is of par- 
ticular interest, a pointing at the Moon (Schmitt 
et al. 1991). The surface brightness of the dark 
side of the moon (the observation was done when 
the Moon was half full) was well above the level 
expected from known sources, and Schmitt et al. 
(1991) suggested that the excess emission was 
bremsstrahlung from solar-wind electrons imping- 
ing on the surface. However, comparing the off- 
Moon intensity of the observation with the cleaned 
RASS data showed a significant excess likely due 
to an LTE. In addition, the enhancement was 
roughly equal to the surface brightness of the dark 
Moon, which suggests a cis-lunar origin for the 
LTEs. 

The study of diffuse X-ray emission from the 
solar system gained interest with the detection of 
X-rays from comets (e.g., Lisse et al. 1996). Af- 
ter some discussion as to the origin of the X-rays 
(e.g., scattered solar X-rays were considered and 
then ruled out) it was eventually determined that 
solar wind charge exchange (SWCX) was respon- 
sible (Cravens 1997). SWCX occurs when a highly 
charged ion of the solar wind interacts with a neu- 
tral atom and picks up an electron in a highly 
excited state. The excited ion then decays to a 
lower energy state emitting a photon with a char- 
acteristic energy of the ion. 

SWCX was subsequently suggested as the 
source mechanism for the LTEs observed in the 
RASS as well as ROSAT pointed observations 
(Cox 1998; Cravens 2000). Because of the time 
scales of their temporal variability measurable by 
ROSAT (a few days or less) the location of the 
LTE production is likely to be relatively local in 
the solar system; either from near Earth’s mag- 
netosheath or within the nearest few AU from 
Earth. 

SWCX has been observed by Chandra (e.g., 
during a dark moon observation, Wargelin et al. 
2004), XMM-Newton (e.g., during an observation 
of the Hubble Deep Field North, HDFN, Snow- 
den, Collier, & Kuntz 2004, hereafter SCK04), 
and Suzaku (e.g., during an observation of the 
north ecliptic pole, Fujimoto et al. 2007). While 
the SWCX emission observed by Chandra clearly 
must originate in the near-Earth environment as 
the pointing was at the dark Moon, that need not 
be the case for the XMM-Newton and Suzaku ob- 
servations. Figure 1 of SCK04 displays the time 


variation of the X-ray emission, solar wind flux, 
and oxygen ionization state ratios for that observa- 
tion while their Figures 2 and 3 show the observed 
spectrum of the SWCX emission. The spectrum 
shows strong lines of O VII and O VIII, as well as 
lines from C VI, Ne IX, and Mg XI, all of which are 
also lines of astrophysical interest. The 0.52-0.75 
keV band light curve has a nearly constant rate 
until a strong solar wind enhancement (as mea- 
sured by the ACE spacecraft) passes the earth, 
and as the solar wind flux and ionization state (as 
measured by the 0 +7 /0 +6 ratio) drop, so does 
the X-ray count rate. On one hand, the nearly 
contemporaneous cutoff suggests a near-Earth ori- 
gin for the SWCX. On the other hand, the nearly 
constant count rate before the drop, in spite of 
the large parameter variations, suggests otherwise. 
Distributing the SWCX emission over a longer 
path length along the line of sight could smooth 
out the light curve, however Collier et al. (2005a) 
demonstrate that the emission would still have to 
originate over a small fraction of an AU. The same 
event has been modeled as due to a large scale he- 
liospheric structure (Koutroumpa et al. 2006), but 
the bulk of the variation is due to local emission. 
On the third hand, the XMM-Newton observation 
geometry was serendipitously good for observing 
SWCX emission from near the sub-solar point of 
Earth’s magnetosheath. This is the region of max- 
imum exospheric emission (Robertson & Cravens 
2003), which again suggests a very local origin. 

Although the event described by SCK04 por- 
trayed the full scale of the problem of the time- 
variable SWCX emission, it was not clear whether 
the emission was due primarily to the solar wind 
interaction with exospheric material in and just 
outside the magnetosheath (a large emissivity over 
a fairly short path-length) or whether the emission 
was primarily due to the solar wind interaction 
with neutral material in the heliosphere (a smaller 
emissivity over a much longer path-length). An 
archival study of sets of repeated XMM-Newton 
observations of the same target showed that time 
variable SWCX can be a problem even when one is 
not observing near the sub-solar point of the mag- 
netosheath, provided that the solar wind flux is 
relatively high (Kuntz & Snowden 2008a) . This is 
made clear by ROSAT and Suzaku observations 
being affected by SWCX as both observatories 
were/ are only able to observe through the flanks 
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of the magnetosheath. (Both satellites had/have 
orbits with altitudes of 550-600 km and were con- 
strained to observe within a range of ~ ±20° from 
the perpendicular to the Earth-Sun line.) How- 
ever, the Kuntz & Snowden (2008a) study did not 
resolve the question of the importance of the mag- 
netospheric origin of the time variable SWCX. 

For the fifth XMM-Newton Announcement of 
Opportunity (AO-5) we proposed for and received 
a 100 ks time critical observation (composed of 
both Guest Observer and calibration time) allow- 
ing us to set the geometry of the pointing to max- 
imize the observed flux from the magnetosheath. 
Despite our best endeavors, and the sterling efforts 
of the XMM-Newton Science Operations Center 
(SOC), we were unable to schedule an enhance- 
ment of the solar wind flux at the time of the ob- 
servation to improve our statistics. However, even 
with a moderate solar wind flux we were able to 
identify a significant correlation between the O VII 
count rate and our model for SWCX emission from 
the near-Earth environment. We were also able to 
use the observation to extract the line intensities 
for O VII and O VIII emission by spectral fitting. 

The goal of this observation was to verify and 
calibrate our model for SWCX emission from 
exospheric material in or outside Earth’s mag- 
netosheath in order to better model and pre- 
dict what, for astrophysical observations, is oc- 
casionally a significant background component. 
This significance was recently made very clear 
by the case of the disappearing warm-hot inter- 
galactic medium (WHIM) near the Coma clus- 
ter. Finoguenov, Briel, & Henry (2003) identi- 
fied excess O VII, O VIII, and Ne IX emission in 
one observation of the many XMM-Newton Coma 
Cluster observations, and this excess emission was 
attributed to the WHIM. However, a recent paper 
by Takei et al. (2008) using Suzaku data from the 
same direction on the sky show no evidence for 
the excess emission. While the solar wind was 
in a relatively quiescent state during the XMM- 
Newton observation we note that the observation 
geometry was similar to the 2001 HDFN observa- 
tion of SCK04 which showed very strong SWCX 
emission. 

Section 2 of this paper describes the X-ray data, 
observation geometry, and ACE solar- wind data, 
§ 3 describes our model for the magnetosheath 
and near-Earth SWCX emission, § 4 describes our 


analysis, and § 5 discusses our results and conclu- 
sions. 

2. Data and Observation Geometry 
2.1. X-ray Data 

The X-ray data used for this study came from 
a single continuous observation obtained with 
XMM-Newton (Jansen et al. 2001; Ehle et al. 
2005) during AO-5, ObsID 0402530201. Table 1 
contains the observation and line-of-sight (LOS) 
details. The observation took place on 2006 June 
4-5 and lasted ~ 96.4 ks. The European Photon 
Imaging Camera (EPIC) detectors (Turner et al. 
2001) were operated in their Full-Frame (MOS1 
and MOS2) and Extended Full Frame (pn) modes 
with medium filters for the MOS detectors and 
the thin filter for the pn. The observation was rel- 
atively unaffected by the soft-proton background 
(see Kuntz & Snowden 2008a). The data were 
processed with XMM-Newton Science Analysis 
Software (SAS 1 ) Version 7.1.2 using the Current 
Calibration Files (CCF) available on 2008 March 
12 . 

The data were reduced using the methods out- 
lined in Snowden et al. (2008), which included 
running the tasks emchain and epchain to create 
calibrated photon event files for the EPIC MOS 
and pn detectors, respectively, espfilt was used 
to filter the data to remove times of soft-proton 
contamination for the spectral analysis. Figure 2 
shows the MOS1 espfilt diagnostic plot with the 
2.5-12.0 keV light curve from the field of view 
(FOV), the 2.5-12.0 keV light curve from the un- 
exposed pixels (pixels in the corners of the de- 
tector that are not exposed to the sky), and the 
histogram of the values in the FOV light curve. 
The diagnostic plots for the MOS2 and pn instru- 
ments are very similar. Spectra were extracted 
from the entire FOV of the MOS data after remov- 
ing point sources detected to a uniform limit of 
2 x 10 -15 ergs cm -2 s -1 and data from CCDs op- 
erating in an anomalous state (Kuntz & Snowden 
2008a). For each point source, the excluded region 
contains 90% of the emission due to the source, 
so the size of the excluded region varies with the 
source’s intensity and distance from the optical 
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axis. The XMM-ESAS 2 software (Snowden & 
Kuntz 2006; Snowden et al. 2008) was used to cre- 
ate model quiescent particle background spectra 
(Kuntz & Snowden 2008a). We did not extract 
a spectrum from the EPIC pn data as we have 
not yet developed a suitable background model for 
that detector. 

We used data from all three instruments for 
the light-curve analysis. We extracted light curves 
in the energy ranges 0.50-0.61 keV (for an O VII 
band) and 2.5-12.0 keV (a soft-proton monitoring 
band) from the MOS1 and MOS2 detectors and 
light curves in the 0.45-0.65 keV and 2.5-12.0 keV 
bands were extracted from the pn (the broader 
O VII band for the pn was required by the slightly 
lower energy resolution). Scatter plots were cre- 
ated of the O VII band versus the hard band (Fig- 
ure 3). For all of the detectors the bulk of the data 
points are strongly clumped, while the outliers are 
linearly correlated. Since the soft proton spectrum 
has a roughly power-law spectral shape affecting 
the entire XMM-Newton EPIC energy range, this 
linear correlation implies that the outliers are due 
to increased soft proton contamination. The count 
rates from all detectors were added together after 
excluding the outliers, i.e., those points not in the 
clumps in Figure 3. Only time periods with ac- 
cepted count rates from all three detectors were 
included in the light-curve analysis. 

2.2. Observation Geometry 

The geometry for this observation had several 
constraints. First, in order to maximize the ob- 
served SWCX emission from the magnetosheath, 
the look direction was required to sweep through 
the subsolar point, which is about 10 R e (Earth 
radii) toward the Sun along the Earth-Sun line 
during typical solar wind conditions. The sweep- 
ing was accomplished by fixing the pointing direc- 
tion for the observation so that the motion of the 
satellite in its orbit would properly position the 
line of sight so it cut through the magnetosheath 
in the desired manner. Because XMM-Newton’s 
apogee distance is also ~ 10 R e and XMM-Newton 
is only able to point within ±20° of the perpen- 
dicular from the Earth-Sun line, we required an 
observation date for which the apogee was near 
the Earth-Sun line (early June). While there was 
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no reason to point at any specific position on the 
sky, in order to maximize the ratio of the SWCX 
emission to the cosmic background we needed to 
choose a direction on the sky where the cosmic 
background is faint in the | keV band. 

After a somewhat painful period of trial and 
error we were able to find an orbit and look 
direction which satisfied the constraints; revolu- 
tion (orbit) #1188, 2006 June 4-5 with a = 
22 /tr 19 m 45.30 s ,5 = +72°26'22.0". The Galac- 
tic column density of H I in this direction is 
~ 3 x 10 21 cm -2 and the RASS | keV count rate is 
~8x 10 -5 counts s -1 arcmin -2 , a very low count 
rate. The RASS \ keV emission is also quite low 
in this direction. Figures 4 and 5 detail the obser- 
vation geometry of the satellite location and look 
direction relative to the magnetosheath. 

2.3. The Solar Wind and ACE Data 

The solar wind, a partially ionized plasma con- 
tinually flowing from the Sun, contains on average 
about 96% protons, 4% alpha particles, and less 
than one percent heavier highly charged species. 
Solar wind densities are typically ~ 5 cm -3 and 
the magnetic field strength is about 7 nT. How- 
ever, solar wind properties are highly time vari- 
able, and over the approximately eleven year so- 
lar cycle the solar wind flux goes from an average 
of about 3.2 x 10 8 cm -2 s _1 at solar minimum to 
about 4.8 x 10 8 cm -2 s _1 at solar maximum, while 
the solar wind speed goes from about 390 km s _1 
at solar minimum to 530 km s _1 at solar maxi- 
mum (Rucinski et al. 1996). The solar wind tends 
to be highly structured with scale lengths as short 
as tens of Earth radii in both the magnetic field 
(Collier et al. 1998) and plasma (Richardson & 
Paularena 2001). At the nominal flow speed of 
about 450 km s _1 , the solar wind takes almost 
four days to flow from the Sun to the Earth and 
over a year to arrive at at the termination shock 
at ~ 100 AU. 

Spacecraft such as the Advanced Composition 
Explorer ( ACE , Stone et al. 1998) and Wind 
(Acuna et al. 1995) monitor the solar wind from 
the LI point, a point about 235 Earth radii up- 
stream where the gravitational force between the 
Earth and Sun balance. Because spacecraft at LI 
typically have halo orbits around this point with 
a radius of about 40 Earth radii, the in-situ ob- 
servations made by upstream monitors may not 
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reflect the exact properties of the solar wind that 
impacts Earth (or another spacecraft) and there 
is a non-negligible lag (~ 1 hour between LI and 
the Earth.) 

The small fraction of heavy high charge-state 
ions in the solar wind produces the soft X-ray 
emission from SWCX. Unfortunately, because of 
their low abundance in the solar wind, measure- 
ments of these minor ion species require special in- 
strumental techniques (e.g., Gloeckler et al. 1995), 
and so are more limited than measurements of the 
proton and alpha particle densities, flow speeds, 
and temperatures. 

During solar minimum, when most XMM- 
Newton observations have occurred, the solar wind 
as observed at Earth has shown a periodic struc- 
ture, oscillating between fast polar coronal hole 
flow and slower interstream flow as the Sun and its 
solar dipole rotate. This behavior is relatively typ- 
ical for observations near solar minimum. At the 
time of the XMM-Newton observation, as shown 
in Figure 6, the solar wind flux measured by the 
ACE spacecraft was usually below its nominal flux 
of 3 x 10 8 cm^ 2 s -1 . The two-hour-averaged 0 + ' 
flux, also from ACE, showed relatively typical val- 
ues. In short, this was a calm solar wind time 
period. 

To give a perspective on the solar wind fluxes 
during the time of the X-ray observation, Figure 7 
shows the integral histograms of the solar wind 
proton and 0 + ' fluxes for the year 2006. The 
range of proton fluxes were typical for the year 
while on average the 0 + ' fluxes were slightly ele- 
vated. 

3. Magnetosheath SWCX Emission Model 

SWCX emission is due to the transfer of an elec- 
tron from a neutral atom, usually H or He, to an 
excited state in a heavy solar wind ion, and the 
subsequent radiative decay. Given the ionization 
states in the solar wind, the resultant photons are 
typically in the extreme UV or soft X-ray region 
of the spectrum. Following Cravens (2000) we de- 
fine the photon production rate (emissivity) for a 
given line as: 

P x = an sw n n (g) (1) 

where n sw is the density of the solar wind protons, 
n n is the density of the neutral “donor” atoms, and 


< g > is their relative velocity. The scale factor a 
contains the abundance of the solar wind ion of in- 
terest, the branching ratios in the radiative decay, 
and the interaction cross-section, which depends 
upon the relative velocity, and the donor species. 
For hydrogen and for medium solar wind activity, 
a reasonable value for an a appropriate for a broad 
X-ray band with E> 100 eV is 6 x 10~ 16 eV cm 2 
(Pepino et al. 2004), and for helium the value is 
about half that of hydrogen. 

The total X-ray emission observed from a spe- 
cific line of sight is given by 

r-S — OO 

4t tI= / P x (s)ds (2) 

Js = 0 

where s is the position along the line of sight. The 
interval over which one needs to integrate is set by 
the phenomena producing the emission. 

For a line of sight originating from a space- 
craft in earth orbit there are two sources of neu- 
tral donor atoms. The first is the extended geo- 
corona, the outermost, most tenuous part of the 
atmosphere, which is composed almost entirely of 
hydrogen. The geocoronal neutral densities are 
obtained from the Hodges (1994) model of the ter- 
restrial exospheric hydrogen density. This model 
is tabulated to 9.7 Re, where the structure is 
fairly smooth; we have extrapolated the model 
at higher R with a R -3 distribution. The sec- 
ond source of donor atoms, which is spatially quite 
complex, is the neutral interstellar medium which 
streams through the heliosphere. The interstel- 
lar medium is comprised primarily of H and He 
where the density of He is about 10% that of H 
in the local interstellar cloud (LIC) which sur- 
rounds the heliosphere (Gloeckler & Geiss 2001; 
Gloeckler et al. 2004). The interstellar medium 
within the heliosphere is initially the same but 
is strongly affected by the Sun. As interstellar 
atomic hydrogen approaches the Sun it experi- 
ences both the attractive force of gravity (oc r~ 2 ) 
and the repulsive force of radiation pressure (also 
oc r -2 ). This creates a hydrogen cavity around 
the Sun with a radius (at the 10% density level) 
of a few AU in the up-interstellar-wind direction 
(A,/? ~ 254° ± 3°, 7° ± 3°) to ~ 10 - 15 AU in 
the downstream direction (Quemerais, Lallement, 
& Bertaux 1993). Because hydrogen is also ion- 
ized as it approaches the Sun, its upwind den- 
sity is significantly greater than its downwind den- 
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sity. Neutral helium is also depleted near the Sun 
due to photoionization, however because the cross 
sections for helium are significantly smaller than 
those for hydrogen, helium can be found much 
closer to the Sun than hydrogen. Further, as 
neutral helium passes the Sun, it is gravitation- 
ally focused, substantially increasing its density 
in the downwind helium-focusing cone. We have 
derived the interstellar neutral densities from the 
Fahr “hot” model (Fahr 1971, 1974). 

There are two important regimes for solar wind 
ion-neutral interactions. The first is the magne- 
tosheath, the region directly behind the Earth’s 
bowshock. The bowshock brakes the solar wind, 
increasing the density of the solar wind in the 
“nose” of the bowshock by roughly a factor of four, 
and by lower factors along the “flanks” of the bow- 
shock. Since the nose of the magnetosheath is also 
the place where the bowshock is the closest to the 
Earth, this is also the magnetosheath region with 
the highest neutral density, and thus the region ex- 
pected to have the greatest SWCX emission. The 
solar wind parameters (density, speed, and tem- 
perature) inside the magnetosheath are given by 
the Spreiter, Sumers, & Alksne (1966) numerical 
global hydrodynamic model. The second regime of 
interest is the free space outside of the bowshock. 
Here there is still a significant density of geocoro- 
nal hydrogen as well as the free-flowing interstellar 
medium. 

The nose of the magnetopause is at ~ 9.7 R# 
under nominal solar wind conditions, and that dis- 
tance varies with the solar wind ram pressure as 
_1 

( n sw u 2 w ) 6 , where u 3W is the solar wind speed. 
Because the exospheric density drops off as R~ 3 
where R is distance from the center of the Earth, 
the exospheric neutral density at the nose of the 
magnetopause varies as ~ R ~ 3 ~ {n sw u 2 sw )~ 2 . 
The soft X-ray emission scales with the product of 
the exospheric neutral density and the solar wind 
flux, n sw u sw , so that the soft X-ray emission goes 

1 3 

as ~ n sw u sw (n sw u 3W )~2 = nj w u 2 aw . The emis- 
sion responds non-linearly to both increases in so- 
lar wind density and solar wind speed varying al- 
most as the square of the solar wind flux. With 
this observation taking place during a relative lull 
in the solar wind, the expected SWCX emission is 
also relatively low. 

We have integrated our model for the interac- 


tion of the solar wind with the geocoronal neutral 
material, including the interactions in the magne- 
tosheath to a distance of 50 Re- We have inte- 
grated our model for the interaction of the solar 
wind with the interstellar material to a distance of 
200 AU; essentially encompassing the entire path 
length through the heliosphere. 

The solar wind conditions enter our model not 
only through the n sw in Equation 1 and through 
the a in Equation 1, which depends on the veloc- 
ity and temperature of the solar wind, but also 
through the Spreiter, Sumers, & Alksne (1966) 
model of the magnetosheath since the size and 
shape of the magnetosheath depends upon the 
strength of the solar wind. We have modeled the 
solar wind as a series of spherical fronts emanat- 
ing from the Sun. The solar wind proton density, 
speed, and temperature for each front was derived 
from the OMNIWeb 3 archive. The data extracted 
from OMNIWeb (King & Papitashvili 2004) had 
a time resolution of one hour, and the data were 
time shifted to the location of the bowshock. The 
proton thermal velocity was derived from the pro- 
ton temperature using 

v t h = \Jzk B T/m p , (3) 

where m p is the proton mass and T is the proton 
temperature, while the total average proton speed 
was calculated using 

(9) = \J Vth, + u lw ( 4 ) 

where u sw is the measured average proton speed. 

As discussed above, the interaction of the so- 
lar wind with the magnetosheath causes mag- 
netosheath to expand and contract. With each 
time step we scale the Spreiter, Sumers, & Alk- 
sne (1966) model of the magnetosheath to the size 
apropriate for the solar wind flux in that time step, 
and then integrate the emissivity along the line 
of sight. Given the limitations of the solar wind 
data, we assume nearly planar propagation near 
the earth and a uniform pressure across the nose of 
the magnetosheath. Such treatment is valid given 
the time-scales involved; due to low count rates 
the X-ray data require a time binning of 15 min- 
utes to achieve a reasonable significance while the 
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time required for a solar wind front to move from 
the nose of the magnetosheath to the level of the 
earth is only about three minutes. Thus, much 
of the dynamic nature of the magnetosheath will 
be averaged out by the time-binning forced by the 
X-ray data. 

4. Analysis 
4. 1 . Light Curve 

In this part of the analysis we consider the 
O VII light curve and test for any correlation 
between the model and the observed count rate. 
The O VII band was chosen as O VII dominates 
the SWCX emission for most observations in the 
XMM-Newton bandpass. The narrow band was 
selected in order to maximize the signal-to-noise 
for the analysis. Note that since the O VII light 
curve was constructed from all of the counts in 
0.5-0.61 keV (MOS) or 0.45-0.65 keV (pn) band, 
the O VII light curve also contains a contribu- 
tion from adjacent spectral bands (due to the fi- 
nite energy resolution of the detectors) as well as 
the underlying spectral continuum. Thus, even 
in the absence of any O VII line emission, the 
O VII band count rate would not be zero. Fig- 
ure 8 shows the O VII band count rate which sums 
the data from all three EPIC instruments. This 
plot also contains scaled curves of the model rates 
for the heliospheric SWCX (dot and dashed line) 
and magnetosheath SWCX (dotted line) emission 
models in order to show their temporal variation. 
(The model emission was arbitrarily scaled for dis- 
play purposes.) The solid curve shows the best- 
fit model which is a constant (to account for the 
non-variable component in the O VII band) plus 
a scaling of the magnetosheath SWCX emission. 
The model heliospheric emission showed too small 
a variation over the observation interval to sig- 
nificantly separate it from a constant component. 
The fit is far from statistically significant with 
y 2 = 191 for 75 degrees of freedom. However, the 
excessive y 2 value appears to be dominated by the 
short-period scatter rather than global deviations 
of the model. If this is the case, the extra scatter 
could possibly be accounted for by additional vari- 
ation in the 0 + ' flux relative to the solar wind pro- 
ton flux which is used in the model. Because the 
heavier ions have much lower abundances the A CE 
solar wind data can achieve reasonable statistics 


only by integrating over longer time intervals. For 
0+' this typically requires time intervals greater 
than an hour. 

Since the y 2 value for the fit is far from ac- 
ceptable we use the correlation coefficient to test 
whether the SWCX model plus a constant can be a 
reasonable representation of the data. The corre- 
lation coefficient is 0.435 for 77 samples, which im- 
plies a probability of < 0.001 that the data could 
be drawn from an uncorrelated distribution (Bev- 
ington 1969). To gain more confidence that the 
correlation is real we considered two additional 
qualitative tests. First, we calculated the y 2 value 
for a range of scale factors. Figure 9 displays the 
resultant y 2 curve with a delta y 2 of ~ 50 between 
the best-fit value and a zero scaling. Thus, our 
magnetospheric model is considerably better than 
a constant model. Similarly, our magnetospheric 
model must be better than the heliospheric model 
alone, which is relatively constant. 

For the second test we examined the variation of 
the y 2 value when the model time is offset from the 
data. Figure 10 shows the result of this test with 
a sharp y 2 minimum which is both deep (a delta 
y 2 value of 30-40) and has a time offset consistent 
with zero. These two tests give us confidence that 
the correlation is real and that the scaling is well 
defined. 

Before fitting the model to the data it was nec- 
essary to scale it to be appropriate for the XMM- 
Newton EPIC detectors for the O VII band. The 
model is in units of keV cm~ 2 s -1 and so must 
be scaled by energy of O VII lines (~ 0.562 keV) 
and the area-solid-angle product for the EPIC de- 
tectors (0.0161 cm~ 2 sr). After this scaling, the 
best-fit scale factor for the magnetosheath emis- 
sion is 1.36, implying that the model underpre- 
dicted the SWCX flux by only ~ 36%. Although 
it is not entirely appropriate, we use the y 2 curve 
shown in Figure 9 to get an idea of the uncer- 
tainty in the magnetosheath model scale factor. 
Using a Ay 2 = 2.71 criteria, which in normal cir- 
cumstances would be 1.6cr, yields an uncertainty 
of 0.32, for a scale factor of 1.4 ± 0.3. 

The average SWCX O VII emission can be ex- 
tracted from the fitted light curve, and is 0.026 
counts s -1 . Assuming that the magnetospheric 
contribution is solely due to the O VII lines, com- 
bining this count rate with the area-solid-angle 
product we find an average intensity of 1.6 ± 
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0.4 photons cm -2 s~ 4 sr -1 (hereafter line units, 
LU), using the same assumption for the uncer- 
tainty as above. 

4.2. Spectral Analysis 

Because of the nature of this observation, long, 
relatively clean, and in a direction with a relatively 
high Galactic column density, it is a good obser- 
vation to study the O VII and O VIII emission of 
the local diffuse X-ray background. 

For the spectral analysis we used Xspec to fit 
the data using a standard model for the emission 
from the diffuse X-ray background. The model 
is comprised of an unabsorbed foreground ther- 
mal component representing the putative Local 
Hot Bubble (LHB) emission, an absorbed thermal 
component, representing emission from the Galac- 
tic disk and halo, and a background power law ab- 
sorbed by the entire Galactic column along the line 
of sight, representing the unresolved extragalactic 
background. The absorption of the thermal com- 
ponent was allowed to vary as the absorption and 
the emission along the LOS may be interleaved. 
The remaining XMM-Newton instrumental back- 
ground was modeled by two lines (A1 K a at 1.49 
keV and Si Ko at 1.75 keV) while the residual 
soft proton contamination was modeled by a power 
law spectrum not folded through the instrumen- 
tal response (Kuntz & Snowden 2008a). In order 
to better constrain the fitted model we used the 
HEASARC X-ray Background Tool 4 to extract a 
RASS spectrum for a circular region with a 1° 
radius around the pointing direction of our obser- 
vation, which was then simultaneously fit with the 
XMM-Newton data. 

APEC 5 models where the relative abundances 
could be individually set were used for the ther- 
mal emission. The abundances for oxygen were 
set to zero and the rest of the abundances were 
fixed at 1. Two additional Gaussians were then 
included to model the O VII emission at 0.56 keV 
and O VIII emission at 0.65 keV. The data were 
then fit to determine the oxygen flux. Figure 11 
shows the data and the best-fit model. To get a 
better perspective on the uncertainties in the spec- 


4 High Energy Astrophysics Science Archive Re- 
search Center, http://xmm.gsfc.nasa.gov/cgi- 

bin/Tools/xraybg/xraybg.pl 

5 http: // cxc. harvard.edu / atomdb / sources_apec. html 


tral fit results for the O VII and O VIII lines we 
used Xspec to create the confidence-contour plot 
shown in Figure 12. 

The fitted values for the fluxes in the O VII and 
O VIII lines are 4.3 ± 0.4 LU and 0.39 ± 0.17 LU, 
respectively. This is a marginal detection of O VIII 
but is significant for O VII. Using the derived value 
for the magnetosheath SWCX O VII from the light 
curve analysis of 1.6 ± 0.4 LU, this leaves 2.7 ± 
0.6 LU originating either as thermal emission from 
the LHB or from the heliospheric SWCX. 

The £ = 111° {£,b ~ 111.14°, 1.11°) line of 
sight described in Kuntz & Snowden (2008b) is 
only 11.8° away from the line of sight described 
here. For that observation, the line of sight passed 
through the flanks of the magnetosheath at an an- 
gle of 106° from the Earth-Sun line. The solar 
wind proton flux was also smaller for that obser- 
vation, 1.39 x 10 8 cm -2 s _1 , and the 0 +7 /0 +6 
abundance ratios were similar for two observa- 
tions. For the £ = 111° line of sight the total 
O VII emission (summed over all emission compo- 
nents) was 3.0^0 4 LU, which is consistent with 
our non-magnetosheath measurement. The to- 
tal O VIII emission in that direction was l.blg' 8 . 
However, it should be noted that the RASS R45 
band value for the £ = 111° LOS is 50 ± 6 x 10 -6 
counts s- 1 arcmin -2 , exceptionally dark due to 
a molecular cloud shadow, while in our LOS it is 
76±7 x 10 6 counts s 1 arcmin 2 . Thus one would 
have expected the i = 111° LOS to be dimmer by 
a third if the underlying spectra were the same in 
both directions; conversely one can attribute the 
difference to SWCX, though the £ = 111° LOS 
should have had a substantially lower magneto- 
spheric contribution. The most likely explanation 
is that the Kuntz & Snowden (2008b) observation 
was unaffected by any significant magnetosheath 
emission and that the excess O VIII line strength 
is due to other thermal emission regions along the 
line of sight in the Galactic disk. 

5. Results and Conclusions 

From our analysis of the X-ray data we derive 
an average O VII flux of 1.6 LU from the magne- 
tosheath and near-Earth SWCX emission. This is 
in very good agreement with the model prediction 
which is 1.2 LU, and provides a very strong valida- 
tion of our SWCX emission model. The rest of the 



observed O VII emission, 2.7 LU, originates fore- 
ground to the dark clouds forming the backstop 
for this observation, and so may originate in the 
heliosphere, LHB, or farther away in the Galactic 
disk. This value agrees reasonably well with the 
3.0 LU results of Kuntz & Snowden (2008b) and 
the 3.5 LU local emission of Smith et al. (2007). 

Our marginal O VIII result of 0.4 LU is sig- 
nificantly lower than those of Kuntz & Snowden 
(2008b) but agrees well with the marginal local re- 
sult of Smith et al. (2007) (0.24 ± 0.10 LU). This 
places an upper limit on O VIII SWCX emission 
from the heliosphere, which should be relatively 
constant over time, of ~ 0.3 LU. 

With a functional model for SWCX emission 
from the magnetosheath and near-Earth environ- 
ment it is now possible to identify those astrophys- 
ical observations which are likely to be affected by 
this “contamination” component. With further re- 
finements of the model and additional calibration 
observations it may be possible to model SWCX 
emission accurately enough to subtract it from the 
data. 

Another interesting possibility suggested by the 
success of this project is the remote monitoring 
of the magnetosheath, magnetopause, and solar 
wind. In the past, observations of the terrestrial 
magnetopause have been mainly restricted to in- 
situ observations. The notable exceptions to this 
are neutral atom imaging (e.g., Taguchi et al. 2004; 
Collier et al. 2005b) and radio plasma imaging 
(e.g., Green & Reinisch 2003; Nagano et al. 2003). 
X-ray imaging, however, may be the best tech- 
nique for global monitoring of the magnetopause 
from, for example, a spacecraft positioned off to 
the side of the Earth in a 1 AU orbit that leads 
or follows the Earth, viewing nearly perpendicu- 
lar to the Earth-Sun line (Robertson & Cravens 
2003; Robertson et al. 2006). Given that all en- 
ergy transfer from the solar wind into the mag- 
netosphere must pass through the magnetopause, 
global imaging of the magnetopause boundary will 
prove invaluable for space weather studies. 
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Fig. 1. — RASS maps of the \ keV dif- 
fuse background before (upper panel) and af- 
ter (lower panel) removal of the LTEs. The 
color bar shows the X-ray intensity in units of 
counts s' 1 arcmin~ 2 . The particle background 
(Snowden et al. 1992; Plucinsky et al. 1993) 
and scattered solar X-ray background (Snowden 
& Freyberg 1993) have been subtracted in both 
maps. 


This 2-column preprint was prepared with the AAS IATgX 
macros v5.2. 
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Fig. 2. — M0S1 diagnostic plot from the SAS 
task espfilt with (upper panel) the histogram of 
the FOV 2.5-12.0 keV light curve, (middle panel) 
the 2.5-12.0 keV light curve from the FOV, and 
(lower panel) the 2.5-12.0 keV light curve from 
the unexposed (to the sky) corners of the detec- 
tor. Data plotted in green indicate accepted time 
periods. In the upper panel the vertical blue lines 
indicate the range over which the histogram data 
were fitted with a Gaussian. The vertical red lines 
indicate the count rate range where the data were 
accepted. 
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2.5-12.0 keV Band Count Rate (counts s _1 ) 



Fig. 3.- O VII band versus hard band count rate 
scatter plots for the MOS1 (upper panel), MOS2 
(middle panel), and pn (bottom panel) detectors. 
The plotted lines show the correlation between the 
outliers indicating that both the hard band and 
O VII band are affected in a similar manner, fur- 
ther indicating that they are due to residual soft- 
proton contamination. 
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GSE-X 

Fig. 4. — Geometry of the observation in the 
GSE X-Z plane, where X points to the Sun and 
Z points toward the north ecliptic pole. The 
curved solid line shows the XMM-Newton orbit 
during the observation while the straight solid 
lines show the look direction. The dotted lines 
indicate the rough location of the inner and outer 
boundaries of Earth’s magnetosheath. The col- 
ored lines provide a rough estimate of the emis- 
sivity ( n sw n n (g )) within the magnetosheath. The 
ragged appearance of the boundary of the mag- 
netosheath emission is due the variability in the 
solar wind strength. 



Fig. 5. — Same as Figure 4 except that the projec- 
tion is in the GSE X-Y plane. The GSE Y axis lies 
in the ecliptic plane and points opposite to Earth’s 
velocity vector, completing the right-handed coor- 
dinate system. 



Fig. 6. ACE solar wind data for the obser- 
vation. The upper curve shows the solar wind 
proton speed in units of km s -1 as measured by 
the SWEPAM instrument (McComas et al. 1998). 
The middle curve is the 0 + ' flux in units of 10 2 
ions cm -2 s -1 as measured by the SWICS instru- 
ment (Gloeckler et al. 1998). The bottom curve 
is the solar wind proton flux in units of 10 8 pro- 
tons cm~ 2 s -1 as measured by the SWEPAM in- 
strument. The vertical lines show the bounds of 
the XMM-Newton observation. 
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Fig. 7. — Integral histogram of selected ACE 
solar wind data for the year 2006. The up- 
per panel (a) shows the solar wind proton flux 
in units of 10 8 protons cm -2 s _1 as measured 
by the SWEPAM instrument. The lower panel 
(b) shows the solar wind 0 + ' flux in units of 
10 5 ions cm^ 2 s -1 as measured by the SWICS in- 
strument. In both panels the vertical lines show 
the bounds of the ranges experienced during the 
XMM-Newton observation. 



Fig. 8.— O VII light curve (data points with 
statistical error bars) and other scaled rates for 
the observation. The solid curve is the best fit 
model which includes a constant plus a scaling of 
the magnetosheath and near-Earth model emis- 
sion (dotted curve) . The dotted and dashed curve 
is a scaled plot of the heliospheric model (SWCX 
with both hydrogen plus helium) emission. The 
time binning for the X-ray data is 15 minutes. 



Magnetosheath Scale Factor 


Fig. 9. — x 2 value as a function of scale factor for 
the magnetosheath SWCX emission model. 
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Fig. 10. — x 2 value as a function of time shift be- 
tween the magnetosheath SWCX emission model 
and data. 
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Fig. 11. — Spectral fit for the XMM-Newton EPIC 
MOS (white and red) and RASS (green) spectra. 
The two lines between 1.5 and 2.0 keV are the 
A1 Ka and Si Ka instrumental lines. The O VII 
line is at 0.56 keV. The O VIII line is not apparent. 



Fig. 12. — x' 2 contour plot. The axis values are 
the Xspec normalizations which are in units of 
counts s -1 arcmin~ 2 . The contours are the 68%, 
90%, and 99% confidence values. 
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Table 1 

Observation details. 


Parameter 

Value 

ObsID 

0402530201 

(a, (5) 

22 hr 19 m 45. 30 s , +72° 26' 22.0" 

(£,b) 

111. 9678°, +12. 8797° 

Start, End Date 

2006-06-04 13:32:04,2006-06-05 16:21:11 

Total, Good Exposure Time 

96.37 ks, 83.34 ks 

Galactic H I 

2.9 x 10 21 cm -2 

RASS R12 

420 x 10 -6 counts s -1 arcmin -2 

RASS R45 

78 x 10 -6 counts s -1 arcmin -2 
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